Is an imaginary fixed point physical or unphysical? 
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It has been proposed that a first-order phase transition driven to happen in the metastable region 
exhibits scaling and universality near an instability point controlled by an instability fixed point of 
a (p 3 theory. However, this fixed point has an imaginary value and the renormalization-group flow 
of the ip 3 coupling diverges at a finite scale. Here combining a momentum-space RG analysis and 
a nucleation theory near the spinodal point, we show that imaginary rather than real values are 
physical counter-intuitively and thus the imaginary fixed point does control the scaling. 

PACS numbers: 05.70.Fh, 64.60.My, 64.60.ae, 64.60.Q- 
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Nowadays, phase transitions are usually classified into 
continuous phase transitions and discontinuous or first- 
order phase transitions (FOPT) [1]. Critical phenomena, 
with the hallmark of scaling and universality [2, 3] exhib- 
ited near a continuous phase transition, have been well 
accounted for by the renormalization-group (RG) theory 
[4, 5]. For a generic FOPT driven beyond its equilibrium 
transition point by an external field in a usual <f> the- 
ory below its critical point, it has been suggested that 
a field-theoretical RG theory of a derived tp 3 model can 
also well account for the universal scaling behavior found 
near such a driven FOPT [6, 7]. A salient feature is that 
here the fixed point of the tp 3 coupling is imaginary in 
value. To reach such an imaginary fixed point, an ini- 
tial coupling with a small imaginary part is sufficient [7] , 
which has been confirmed in a nonpcrturbative RG anal- 
ysis [8]. However, for a physical coupling of purely real 
initial value, the RG flow has to diverge at a finite scale 
in order to achieve an imaginary part. This leads to con- 
cern as to whether the perturbation analysis is still valid 
and thus whether such a theory can indeed describe the 
scaling behavior of the FOPT or not. 

To answer these questions, we need to know what phys- 
ically this finite scale is, why it leads to divergence, and 
where the imaginary part comes from. Here, combining 
a momentum- shell integration RG analysis and a nucle- 
ation theory [9] near the spinodal point [10, 11], we show 
that, at the finite scale, the free-energy cost for nucleation 
out of the metastable state in which the system lies and 
the metastable well itself vanish. This places the system 
exactly at a true instability point and exhibits the diver- 
gence. The integration for the partition function then 
diverges and has to be analytically continued to complex 
plane in order to be physically meaningful. As a result, 
the system enters the imaginary domain and can thus 
reach the imaginary fixed point. Therefore, counter the 
intuition that only real values are physical, here imagi- 
nary values are physical instead. The divergence at the 
characteristic scale just signals the irrelevant degrees of 
freedom have been eliminated and the true instability 
point has been reached and thus does not spoil the per- 



turbation RG analysis. 

First we briefly derive the tp 3 model. Consider the 
usual 4 model with a Ginzburg-Landau Hamiltonian 



Hi[4>] = / dx 
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rP + -rM + ^K{Vcj>Y - H<t> 
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for a scalar order parameter <f> in the presence of its con- 
jugate field H , where r is the reduced temperature and 
g a coupling constant. The constant K has also been 
introduced for the later RG study. Although metasta- 
bility is a dynamic process, statics is sufficient for our 
present purpose. In the mean-field approximation, it is 
well known that H4 has a critical point at r = and 
H = 0. On the other hand, for each r < 0, there is an 
FOPT between two equilibrium phases, each with an or- 
der parameter of opposite sign to the other. However, 
this transition can only take place beyond the spinodal 
point at which the barrier between the two phases van- 
ishes in the same approximation. In order to study the 
property of the transition, we shift the order parameter 
by <f> — Mo + <P with a uniform Mo, neglect the irrelevant 
terms, and arrive at the (p 3 model [6, 7] 



H[<p] 



dx 



1 



1 



(2) 



which has the spinodal at r = r + = and 

h = H — rMo — <?Mo/3! = as its mean-field insta- 
bility point, where v = gM . Given the partition 
function is obtained by a functional integration over all 
configurations 



Z = J Vipcxp(-rl) 



(3) 



and the free energy is F = — In Z. In the following, we 
shall concentrate on the model at h = 0. In fact, a shift 
of (p can change h to r and vice versa and only one of 
them is independent. 

In the mean-field approximation in which tp assumes a 
uniform most probable configuration M, the free energy 
density is 



f(M) - F/V 
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gM 3 



(4) 
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FIG. 1. (Color online) (a) <^ 3 mean-field free-energy density 
/ for different r. (b) The potential V = — / for r < 0. 



where V is the system volume. As seen in Fig. 1(a), for 
t ^ 0, there is a metastable minimum and a saddle point 
at M m = —2t/v and M s = 0, respectively, for r < 
with a difference in free-energy density A/ = 2|t| 3 /3i> 2 . 
When r > similar behavior appears; only their roles 
exchange. So, we shall assume r < for definiteness as 
we shall see below that its fixed point value is of this sign. 
The transition from M m and M s happens at r = 0, the 
spinodal point, at which both states merge. This corre- 
sponds to the real system moving from its own metastable 
to the spinodal point, at which the real FOPT takes place 
corresponding in turn to the transition from M s to the 
unbounded left side in Fig. 1(a). Although the two tran- 
sitions appear different, they occur at the same spinodal 
and it can be shown that in its vicinity the mean-field 
exponents describe correctly the real FOPT [7]. This 
supports strongly the relevancy of the (p 3 theory. 

Next we present the momentum-space RG study of 
the model, which also serves as an introduction to the 
problem. Comparing T~L with H 4 , one sees that the role 
played by the spinodal point to the ip 3 model is similar 
to that played by the critical point to the 4> 4 theory. This 
prompts an RG study to take into account fluctuations. 
The usual RG consists in first integrating out the order 
parameter with wave numbers between the cutoff A and 
A/s with a scaling factor s > 1 and then rescaling the 
length to restore the cutoff. The remaining order pa- 
rameter has also to be rescaled accordingly in order to 
ensure invariant [4, 5]. A previous RG study of a model 
with both (f> 4 and (j> 3 couplings has been performed to 
investigate pseudocritical phenomena near the spinodal 
point, with a conclusion that a true spinodal does not ex- 
ist [12]. Collecting those terms relevant to H and setting 
s = 1 + SI, we find the differential RG equations [12] 



-m = - r > K 
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(5a) 
(5b) 
(5c) 



where 77 is an exponent, d the space dimensionality, 
e = 6 - d, and A n = K d K d - 6 K n / 2 with K d = 



2/[(47r) d / 2 r(d/2)] (r is the gamma function). The flow 
equation (5) has, besides the Gaussian fixed point at 
t = and v — 0, a nontrivial fixed point 

2d A 
4-d 



,*2 



(12 + d)A 6 
i-d 

'16 



2d 



A 6 v* 2 



12 + d 



-A 4 v* 2 = - 



AA 4 



(6a) 
(6b) 
(6c) 

,*2 — 



which is, indeed, the "instability" fixed point with u* 
v* 2 Aq = — 2e/3 and the "instability" exponent rj = —e/9 
to one-loop order found by a field-theoretical RG method 
[6, 7] when d — 6, the dimension on which the e expansion 
is based. Although the fixed point is imaginary in value, 
one can show from Eqs. (5) and (6) that it is infrared 
stable for d < 6 similar to the fixed point for critical 
phenomena and thus controls the large-scale behavior. 
There exist also evidences that this fixed point influences 
the flows of effective exponents in a Monte Carlo RG 
study of the FOPTs in the two-dimensional Potts models 
[13]. The reason that we call it instability fixed point will 
become clear below. 

In order to see how to reach the imaginary fixed point, 
we now solve the flow equation (5). As AT remains a 
constant, Eq. (5b) leads to 



v 2 {l) 



v 2 Ae 



(7) 



(Ae + A 6 vl)e- d - A 6 v 2 ' 

where vq is the initial coupling constant at I = as can 
be easily checked. For e < or d > 6, This solution 
approaches the Gaussian fixed point continuously as I — » 
00. For e > 0, v 2 also approaches v* 2 correctly in the 
same limit. However, the course of the approach varies. 
Whereas for a complex vq , the fixed point is continuously 



reached, for a real Vq, 
satisfying exp(— el s ) = 

L 



v diverges at a characteristic scale 
: A 6 v 2 /(Ae + A 6 v 2 ), or, to O(l), 

A AK ' , (S) 



K d A d - 



AqVq -fVdil Vq 

where we have kept A for dimensional reason, since l a is 
dimensionlcss. Larger than this scale, v becomes purely 
imaginary and can thus reach the fixed point. From 
Eq. (7), the solution to Eq. (5c) can be expressed as 



r(i) = r e 2/ + 



for I < L and 
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for / > l s , where To = t(1 = 0) and $(z,fc, a) = 
S^Lo z ™/( ri + a ) k 1Si ^ ne Hurwitz-Lerch transcendent, 
which has singularity when a equals negative integers 
and branch cut discontinuity in the complex z plane run- 
ning from +1 to oo. So t diverges at l s too. Because 
$(0, 1, a) = 1/a, t(1) approaches its fixed point value (6c) 
correctly as I — > oo. However, to reach this, To has to as- 
sume a particular value that nullifies the terms in the 
braces in Eq. (9b), i.e., to cancel roughly the accrual up 
to l s . This may appear strange as to may be supposed 
to be zero to suppress the growth of t as it is relevant. 
However, it can be shown that in the 4> 4 model for criti- 
cal phenomena, t has also to take a particular value to 
cancel the value at I = of a function similar to the last 
one in Eq. (9b). 

Our primary task is then to show what the charac- 
teristic scale is. To this end, we employ the theory of 
nucleation near the spinodal point [10, 11], which is the 
extension of the nucleation theory near the coexistence 
curve [9] to the spinodal point in the case of long-range 
interactions. The only difference is that we consider t 
instead of h studied there. According to these theory 
[9-11, 14], the critically nucleating profile is given by the 
saddle-point solution of the Eulcr-Lagrange equation 



i 



(10) 



for H. Upon assuming that the profile depends only on 
the radius denoted by x and neglecting the first derivative 
[10, 11], Eq. (10) becomes 



d 2 (p s 
dx 2 



2 _ 



dV 

dip s 



(11) 



with a potential V = ~f(ip s ). Equation (11) describes 
a particle of unit mass moving in the potential V with 
tp s and x representing its displacement and time, re- 
spectively. Therefore, there is a bounded solution if the 
particle has an initial total energy equal to V(t/v) — 
V(—2t/v) as shown in Fig. 1(b). The solution is 



2t 

<Ps = 

v 




(12) 



which can be checked to have correct limits. This intro- 
duces a correlation length £ ~ J K/\r\ that diverges as 
t — )• 0. However, it does not relate to l s - Note that the 
neglected first derivative acts as a frictional force and has 
been shown to change only the shape of the profile [11]. 
In particular, a bound state exists apparently in this case 
too and the form of the correlation length should keep up 
to a constant. To create the critical nucleus, there is a 



free-energy cost [9-11, 14] 
AF = F{ Vs ) - F(M m ) 
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where uses have been made of Eqs. (11) and (12) and 
we have neglected higher-order corrections to ip s . Note 
that AF ~ £ d A/ reasonably and thus its form as given 
in Eq. (13) must be valid even the neglected first deriva- 
tive is taken into account. Equation (13) shows that for 
d > 6, the free-energy cost for nucleation diverges as the 
spinodal point is approached; while for d < 6, AF van- 
ishes with t. This has been argued against existence of 
the spinodal [10, 11] but serves as a physical origin of the 
borderline critical dimension for the e expansion. 

Now it can be seen that AF is related to the charac- 
teristic scale l s - Apparently, one sees that in d = 6 on 
which the e expansion is based, the K 3 /v 2 dependence 
of the free-energy cost AF just coincides with l s . Note 
that except for possible dimensionless multipliers it is in 
fact the only combination that can have the dimension 
of A d ~ 6 . Indeed, if one assumes that the dimension of ip 
is £, denoted as [ip] = £, i.e., <p ~ A^, then [r] = d — 2Q, 
[v] = d — 3£, and [K] = d — 2 — 2£, because H and F are 
dimensionless as we have absorbed the thermal energy 
into their definitions. So, [F 3 /w 2 ] = d — 6. 

More explicitly, we can show formally that l s is just 
the scale at which AF vanishes. Indeed, using Eq. (5), 
we find 
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where we have kept e on the left hand side for definiteness 
but expanded it on the right hand side. Equation (14) 
indicates that W 2 /i; 2 reduces constantly upon coarse 
graining. Accordingly, we have 
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(15) 



Therefore, to O(l), exactly at l s given by Eq. (8), 
t c / 2 /Aqv 2 vanishes. One can check that a constant mul- 
tiplier does not of course change the result and thus AF 
indeed vanishes at l s . 

Equation (15) shows, in fact, the effect of v~ 2 and 
hence the effect of the divergence. Note that Eq. (14) 
is only formal since the neglected terms diverge at l s . 
However, its right hand side only contains t, while the 
flow equation for t £ / 2 itself, for example, contains, be- 
sides those ts in Eq. (14), also v 2 , which is the source 
of the divergence. Also, the flow equation for r/v is not 
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simple either. This seems to indicate that Eq. (14) is 
reasonable. Of course, if we ignore simply -W 2 with the 
e expansion, Eq. (15) is just the leading-order terms of 
v~ 2 as can be checked from Eq. (7). So, at least to the 
leading order in the e expansion, AF indeed vanishes at 
l s . In fact, since v~ 2 equals zero exactly at l s , all other 
quantities multiplying v~ k with k > vanish at l s even 
though they themselves diverge there. Accordingly, at 
l s , not only AF, but also t/v vanish as can be checked 
numerically from the solutions (7) and (9). Therefore, at 
I = l s , the shape of the unform part of the Hamiltonian 
becomes similar to the chain line with r = in Fig. 1(a). 
Moreover, since v diverges, it becomes extremely steep. 
In other words, due to the divergence, the system be- 
comes extremely unstable at Z s . 

Knowing what l s and the effects of the divergence are, 
we can then understand their related physics and why 
the flow becomes imaginary. Given a i>o, at least near 
d = 6, the l s divides all the modes into two sets. Whereas 
the short-wavelength modes that have their wavelengths 
shorter than that which corresponds to l s have a finite 
frcc-cncrgy cost AF for nucleation, the long-wavelength 
ones have none. As the RG e expansion is about the 
mean-field theory in which no fluctuations exist, a finite 
AF thus prevents such modes from escaping the free- 
energy well. Accordingly, the short-wavelength modes 
probe basically the local metastable state in Fig. 1(a) 
and do not know the existence of the other stable, al- 
beit unbounded, state due to their short-range fluctua- 
tions. As a result, they cannot trigger the FOPT and 
are thus irrelevant to it. On the other hand, the other 
modes are extremely unstable as they lie just on the true 
instability point and are clearly relevant to the FOPT. 
According to the theories of nucleation [9-11], the free 
energy of a metastable state becomes complex when the 
saddle point contributes, because the latter is unstable to 
perturbation towards the stable state and thus analytical 
continuation is needed in obtaining Z. Thus, although 
the short- wavelength degrees of freedom can be safely in- 
tegrated away, the contour for integrating the unstable 
long- wavelength ones has to be deformed analytically to 
obtain a meaningful convergent free energy. Indeed, such 
an integration in Eq. (3) diverges for ip < 0. It has ac- 
cordingly to be performed along a contour such that the 
real part of ip 3 is positive. This can be done by integrat- 
ing along the imaginary axis for ip < 0, which amounts 
to replace ip with iip and thus v becomes iv along the 
new contour as speculated [7]. Therefore, the RG flow 
of real values corresponds to elimination of the irrelevant 
degrees of freedom in order to place the system at the 
true instability point, which exhibits divergence; while 
the flow beyond becomes imaginary in order to be physi- 
cal counter-intuitively and drives the system to the fixed 
point. Note that as the division of modes is set by AF, 
we may well reverse the above derivation: It is physically 
more plausible that the vanishing of AF at l s results in 



a divergent v 2 from Eq. (15). 

The above idea can be put in an explicit way. Taking 
into account the convergence of the ip 3 model and the 
division of different degrees of freedom, we may write 
from the beginning the ip 3 interaction of H. as 




with two real parameters i>i and V2, where the two in- 
tegrations are over spatial extents less and greater, re- 
spectively, than that which corresponds to the scale l s . 
In this paramctrization, the first term is irrelevant and 
disappears beyond l s upon coarse graining and can thus 
be ignored; while the second one is the only relevant in- 
teraction for a meaningful convergent description of the 
real FOPT from the instability point to the left side 
in Fig. 1(a). The relevant model coincides then with 
the one proposed for the Yang-Lee edge singularity [15]. 
Its RG behaves benign as it has no divergences and a 
real fixed point that is just the imaginary one found 
above. This therefore confirms the perturbation anal- 
ysis and the physics of the imaginary fixed point . More- 
over, the model justifies directly without the mapping 
employed previously [7] that the edge singularity and the 
FOPTs fall in the same universality class as they now 
share the same relevant model. However, the real Hamil- 
tonian (2) itself can automatically find its proper fixed 
point through divergence. Note that it appears that nu- 
cleation may also be governed by the same fixed point. 

Concluding, we have shown that, at the characteristic 
scale at which the RG flow diverges, the free-energy cost 
for nucleation out of the metastable state vanishes and 
the system becomes extremely unstable and then enters 
an imaginary domain for convergence. One may still seek 
how the perturbation expansion in d = 6 in which AF 
does not depend on r can be carried over to other di- 
mensions. However, the clear physical picture emerges 
should suffice to show the imaginary fixed point of the 
tp 3 model does be reachable physically and control the 
scaling and universality of driven FOPTs. 

I thank Shuai Yin, Yantao Li, Guangyao Li, and Ning 
Liang for their useful discussions. This work was sup- 
ported by NNSFC (Grant No. 10625420) and FRFCUC. 



[1] M. E. Fisher, Rep. Prog. Phys. 30, 615 (1967). 

[2] H. E. Stanley, Rev. Mod. Phys. 71, S358 (1999). 

[3] For a recent review, see M. Barmatz and I. Hahn, Rev. 
Mod. Phys. 79, 1 (2007). 

[4] K. G. Wilson and J. Kogut, Phys. Rep. C 12, 75 (1974). 

[5] S. K. Ma, Modern Theory of Critical Phenomena (Ben- 
jamin, Reading, MA, 1976). 

[6] F. Zhong and Q. Z. Chen, Phys. Rev. Lett.95, 175701 
(2005). 

[7] F. Zhong, arXivl205.1400. 



5 



[8] Y. Li and F. Zhong, arXivllll.1573. 

[9] J. S. Langer, Ann. Phys. (N. Y.) 41, 108 (1967). 

[10] W. Klein and C. Unger, Phys. Rev. B28, 445 (1983). 

[11] C. Unger and W. Klein, Phys. Rev. B29, 2698 (1984). 



[12] Y. Saito, Prog. Theor. Phys. 59, 375 (1978). 

[13] S. Fan and F. Zhong, J. Stat. Phys. 143, 1136 (2011). 

[14] J. W. Cahn and J. E. Hilliard, J. Chem. Phys. 28, 258 

(1958); 31, 688 (1959). 
[15] M. E. Fisher, Phys. Rev. Lett.40, 1610 (1978). 



